Global Deep-MOND Parameter as a Theory Discriminant 
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Different formulations of MOND predict somewhat different rotation curves for the same mass 
distribution. Here I consider a global attribute of the rotation curve that might provide a convenient 
discriminant between theories when applied to isolated, pure-disk galaxies that are everywhere deep 
in the MOND regime. This parameter is Q = (V 2 )/V^, where (V 2 ) = M" 1 / 2nrE(r)V 2 (r)dr, with 
E(r) the disk's surface density, M its total mass, and Voo is the asymptotic (constant) rotational 
speed. The comparison between the observed and predicted values of Q is oblivious to the distance, 
the inclination, the mass, and the size of the disk, and to the form of the interpolating function. 
For the known modified-gravity theories Q is predicted to be a universal constant [independent of 
E(r)]: Q = 2/3. The predicted Q value for modified-inertia theories does depend on the form of 
E. However, surprisingly, I find here that it varies only little among a very wide range of mass 
distributions, Q « 0.73 ± 0.01. While the difference between the theories amounts to only about 
5% in the predicted rms velocity, a good enough sample of galaxies may provide the first discerning 
test between the two classes of theories. 

PACS numbers: 04.50.-h 98.52.Eh 98.80.-k 



I. INTRODUCTION 

In applying modified Newtonian dynamics (MOND), 
we use several formulations and theories (for a recent re- 
view of MOND and its formulations see [1]). The nonrel- 
ativistic formulations in use are the nonlinear version of 
the Poisson equation [2], quasilinear MOND (QUMOND) 
[3]-both of which may be classified as modified-gravity 
(MG) theories-and a class of so-called 'modified-inertia' 
(MI) theories [4, 5]. 

We do not know yet which of these formulations, if 
any, is in the right direction. It is important to try and 
decide between these options, because they each point 
to a different direction in constructing relativistic formu- 
lations, and, more generally, in pinpointing the deeper 
origins of MOND. Of the known relativistic formulations 
of MOND, TeVeS [6], MOND adaptations of Einstein- 
Aether theories [7], most versions of bimetric MOND 
(BIMOND) [8], and those based on a polarizablc medium 
[9] have a version of the nonlinear Poisson theory as their 
nonrelativistic (NR) limit; while certain versions of BI- 
MOND tend to QUMOND in the NR limit. 

But, to my knowledge, no systematic endeavor has 
been undertaken to observationally discriminate between 
the different theories. The main reason for this must be 
that despite the pronounced conceptual differences be- 
tween them, the differences between their salient predic- 
tions are small: by and large, these predictions follow 
directly from the basic tenets of MOND, which all for- 
mulations share. These tenets are [10]: departure from 
Newtonian dynamics for accelerations a < ao, and space- 
time scale invariance in the deep-MOND limit (DML) 
aCao, where ao is the MOND acceleration constant. 

Arguably, the most effective way to distinguish obser- 
vationally between the theories would be based on their 
somewhat different predictions of rotation curves (RCs) 
of disk galaxies. This, however, is not so easy. For all 
galaxies, the basic tenets alone dictate that the rotational 



speed becomes constant at large radii, and they also fix 
the value of this speed. For galaxies with high acceler- 
ations near their centers, all theories predict Newtonian 
speeds there, and the basic tenets predict that the tran- 
sition from Newtonian to MOND behavior should occur 
around the radius where V 2 (r)/r = ay. So the predicted 
curves follow almost in full from only the basic tenets. 
For galaxies with a <C ao everywhere, we shall see below 
that the differences between theories of the predicted rms 
values of the rotational speed are only about 5 percent. 
Indeed, Ref. [11] showed through a few numerical exam- 
ples that the predicted differences between the nonlinear- 
Poisson formulation of MG, and MI theories, are of this 
order. 

In addition, some of the differences between the predic- 
tions may be reduced by choosing different forms of the 
interpolating function that appears in all theories, and 
also by assuming different stellar mass-to-light (M/L) 
values (which are not known accurately) when predicting 
the RC of a given galaxy. 

Because of these hindrances, one would require a large 
sample of galaxies with very good quality data to conduct 
meaningful tests. But the task of computing RCs for 
many galaxies in MG theories is quite onerous. Perhaps 
because of all the above obstacles, no attempt has been 
made so far to conduct a systematic discriminating study. 
All but a few MOND RC analyses of real galaxies to 
date were performed using the rather more manageable 
MI predictions. An exception is the recent analysis of 
several RCs using QUMOND [12]. 

Here, I discuss a useful shortcut for performing such 
an analysis using not the full RC, but a global parameter 
predicted to attain different values by different theories. 
Its main advantage is that it eliminated the need to com- 
pute RCs for different galaxies, since, as I show here, it 
takes an almost universal value within each theory class. 
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II. THE Q PARAMETER 

It follows from the basic tenets of MOND, and hence in 
any present MOND theory, that the asymptotic circular 
rotation speed around any isolated mass is constant and 
depends only on the total mass M: 

V x = {MGa ) 1/4 . (1) 

(Strictly speaking, only a proportionality is dictated, but 
a is normalized so as to give equality.) Consider a pure- 
disk galaxy of surface density £(r) and RC V(r), that is 
wholly in the deep-MOND limit (DML): V(r) 2 /r < a 
everywhere. We define 

Q ^^l = ^Vl[ 2 ^{r)V 2 {r)dr (2) 

({V 2 ) here is not to be confused with the velocity dis- 
persion). When testing MOND itself it behooves us to 
work with the velocities themselves since the amplitude 
of the RC is also a prediction of MOND. But since all 
MOND theories predict the same asymptotic speed, we 
gain much, when only comparing MOND theories, by 
considering the normalized speeds as in the definition 
of Q. I restrict myself to pure disks because otherwise 
the unknown kinematics of the bulge component would 
enter. Restriction to the DML has the following advan- 
tages: (i) It frees us from dependence on the exact form 
of the interpolating function, since the DML form of this 
function is the same for all theories, (ii) Because of the 
scaling properties of the DML, and the choice of normal- 
ized speeds in Q, its predicted values depend only on the 
form of £(r), not on the mass and size of the galaxy, pre- 
dicting the same value for all surface densities nT,(r/h) 
independent of k and h. This means that comparison of 
the predictions with the data is oblivious to the assumed 
distance to the galaxy or its (overall) inclination. This is 
not the case for disks that are not in the DML. (iii) In 
this limit the known MG formulations of MOND predict 
a universal value of Q, saving us the need for onerous nu- 
merical solutions of the nonlinear potential problem for 
each galaxy separately, (iv) A surprising result of the 
present work is that the Q value predicted for MI the- 
ories is also almost universal (while different from that 
predicted by MG theories). 

A. The predicted Q value— modified-gravity theories 

The two nonrelativistic, MG theories in use today, the 
nonlinear Poisson formulation and QUMOND probably 
predict somewhat different RCs for the same mass dis- 
tribution (although I am not aware of any actual com- 
parison between the two). However, they predict the 
same DML virial relation: for any isolated, stationary, 
self- gravitating DML system made of masses m i} both 
theories predict [3, 13] 

<V 2 ) = |(MG ao ) 1/2 (l-£^ 3/2 ), (3) 



where (V 2 ) = M^^^^niiV 2 is the mass weighted mean- 
squared (three-dimensional) velocity in the system (the 
velocities are measured with respect to the center of 
mass), and qi = nii/M are ratios of particle masses to 
the total mass. 1,2 Applying this to a disk galaxy, where 
mi are stars or gas clouds, we have J2<k^ 2 ^ 1 ; if tnc 
system is made of N masses with q^ ~ 1/N then q'J 2 ~ 
N' 1 / 2 . We thus have from Eqs. (l)-(3), in the large-iV 
limit, 

Q = 2/3, (4) 

with a correction of order A -1 / 2 . 

So both theories predict the same rms velocity (and the 
same asymptotic speed) for DML galaxies, which must 
mean that their predicted RCs are very similar. 

B. The predicted Q value— modifled-inertia theories 

We do not yet have an example of a full-fledged MI 
theory [4, 5], but, interestingly, we do know exactly the 
form of the predicted RC for all such theories. A general 
theorem [4] states that for circular orbits in axisymmet- 
ric potentials (applicable to the circular motions in disk 
galaxies) the MOND acceleration, g(r), and the Newto- 
nian one, g N (r), at radius r, are related algebraically: 
g(r) = g N (r)v[g N (r)/ao], where v(y) is a function that 
depends only on the theory (it is derived from the action 
of the theory restricted to circular orbits). The basic 
tenets of MOND dictate the asymptotic behavior of v(y) 
as small and large arguments. In particular, they dictate 
the unique DML relation universal to all MI theories: 

g(r) = [oo^r)] 1 ^. (5) 

(For some disks g N may point outward at some r; then 
g does as well; but I ignore this for simplicity.) So, the 
predicted value of Q can be written in terms of the New- 
tonian RC of the disk, V/v(r): 

o roo 

Q -Wm^L ^MWr)*. (6) 

This does not lead to a general virial relation of the type 
of Eq.(4): the value of Q is not universal. 

To sec if Q can serve as a useful discriminant, I cal- 
culated its value for various classes of S(r) to see how 



1 In calculating (V 2 ) the internal motions within the masses rrn 
are immaterial. If we want to include those we need to consider 
the constituents of ra; as elementary masses, and ensure that 
their internal dynamics are in the DML. If they are not (as in 
stars), or if m; themselves are not even sclf-gravitating (as in 
atoms), we need to consider them as elementary bodies. 

2 It does not matter if the accelerations inside or in the near vicin- 
ity of mi are high. It is only required that the 'mean-field' ac- 
celerations are small everywhere. 
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similar they are to each other, and how different they are 
from the universal Q = 2/3 of MG theories. 

For Kuzmin disks having £(r) oc [l + (r//i) 2 ]~ 3 / 2 I find 



Q = (16/15)(2/^) 1 / 2 [r(7/4)] 2 » 0.719. 



(7) 



A family of finite galactic disks has been described 
in Ref. [14] having a Newtonian RC Vjv(r) = r\u a 
{u = r/r a , rj = [n 1 / 2 (2a + l)T(a + l/2)/2T(a + 
l)] 1//2 (A'/G/r ) 1 / 2 } within the material disk, whose sur- 
face density vanishes at and beyond r , and inside r : 

(2a + l)M , , n 1 3 „ 

= 27 rr 2 ^ (1 ~ V ] 2Fl{l - a >2 ] 2 ; V } ' 

(8) 

where u = (1 — u ) ' , and 2-F1 ar e hypergeometric func- 
tions. For these, I find 



7r 3 / 2 (2a+l)r(q + l/2) 
18 r(a+l) 



,1/2 



r(a/2 + 5/4) 
r(a/2 + 7/4) 



• (9) 



As a varies in the relevant range 3 1/2 < a < 3/2, Q 
varies only between 0.740 and 0.726. Some interesting 
special cases are the 'isodynamic' disk [15], which has a 
constant acceleration (both in Newtonian and MI DML), 
with a = 1/2, for which Q = n/2 1 / 2 3 w 0.740; the 
Kalnajs disk, with a = 1, for which Q w 0.736; and 
a = 3/2, corresponding to rigid rotation in MOND, for 
which Q = 2 7 / 2 /3 5 / 2 w 0.726. 

For all exponential disks, which have £(r) oc e~ r / h , 
one has 



- 8 jT 

Jo 



S 5 / 2 e- 2s [/ (s)^o(s)-/i( S )^i( S )] 1/2 ^, (10) 



where and K v are the modified Bcssel functions of the 
first and second kind (s = r/2h). Evaluated numerically, 
Q w 0.733. 

To expand the range of disk models, I also considered 
double exponentials with 



£(r) cx e- r/h + pe~ qr ' h , 



(11) 



with the two disks having a mass ratio of (3/q 2 . For these, 
Q can be written as 



Q 



q 2 + PJo 



^ 2 (e- 2s + Pe- 2 « s ){I ( S )K { S ) 



-h{s)Ki{8) + /3q[I (qs)K {qs) - I^K^qs^^ds. 

(12) 

For the pairs (q,/3) = (2,1), (5,1), (2,4), (3,9) (the 
last two correspond to equal masses), I find Q ~ 



0.7312, 0.7319, 0.7299, 0.7257. Taking -1 < /3 < 0, 
we get exponential disks suppressed near the center 
(an actual hole occurs if q\/3\ > 1, but if q\(3\ is too 
large becomes negative). For the pairs (q,/3) = 
(2,-0.5), (1.5,-0.7), (1.2,-0.8), (2,-0.55), (4,-0.3), 
Q w 0.7340, 0.7346, 0.7346, 0.7341, 0.7332. 

Surprisingly, and for reasons that I do not understand, 
all the above disk models (which are all that I tried) give 
very near values of Q = 0.73 ±0.01, compared with the 
universal MG value of Q = 2/3. The difference is small, 
amounting to about a 5% difference in the rms veloc- 
ity, indicating that the predicted RCs in all the theories 
concerned differ only a little from each other. 4 Using 
Q as discriminator may be rather demanding, consider- 
ing the possible sources of systematic errors (see section 
III). Still, hopefully, with accurate enough data for a 
large enough sample of galaxies, these differences might 
be used to distinguish between the two classes of theories. 



III. DISCUSSION 

Since the differences in rms velocities predicted by the 
two theory classes are only about 5%, the proposed test 
will be potentially stymied by measurement and various 
systematic errors. The main concerns I can think of are 
as follows: (i) Accurate determination of for low- 
acceleration galaxies is not always possible, since their 
speed at the furthest observed radius is, typically, larger 
than by a few percent. For example, the predicted 
DML RC of an exponential disk peaks at about 5 scale 
lengths, beyond which it declines to V x , which is about 
5% lower. This fact has been discussed in Ref. [16] in 
connection with its effect on the measured intercept of 
the M — V x relation [Eq.(l)]. (ii) There may be a possi- 
ble unaccounted for contribution to the integral in Q from 
radii beyond the furthest observed one. Since V ~ 
there, the missing contribution to Q is the fractional disk 
mass beyond this radius, (iii) In very low velocity galax- 
ies, which some low-acceleration disks are, the analysis is 
susceptible to uncertain, asymmetric-drift corrections. It 
is thus advisable to exclude such galaxies from the anal- 
ysis, (iv) In determining Q observationally, we need the 
observed RC, as well as £(r) (both only up to an imma- 
terial normalization). If the disk is made of gas only, or 
of stars only (for which we may assume a constant M / L 
ratio) this can be cleanly done. But if both components 
contribute, the exact shape of E(r) depends on the stel- 
lar M/L value, which is not known a priori with enough 
accuracy. It is thus advisable to use galaxies that are 
dominated by either component (as used in Ref. [17] for 
a clean test of the MOND relation Eq. (1)). 



3 Disks with a < 1/2, such as the Mestel disk (with a = 0, for 
which Q 0.730) are not legitimate DML disks, because their 
accelerations diverge at the center. For a > 3/2 the MOND RC 
is (unrealistically) concave. 



4 Looking at the few comparisons of full RCs in Ref. [11] we sec 
that the contribution to the difference in Q comes from the bulk 
of the disk (roughly around where S takes half its central value). 
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Other possible systcmatics to worry about are, depar- 
tures from the DML , and the applicability of the razor- 
thin-disk approximation, assumed all along. 

We may also consider other moments of the velocity 
curves as discriminants. While they may have some ad- 
vantages, they will have the drawback that there is no 
simple expression for them (not even as integrals) in MG 
theories, and they will have to be computed for each sur- 
face density law separately (for MI theories they can be 



calculated with the same ease as Q). For example, I 
looked at Q 4 = (V 4 ) 1 ' 2 . For Kuzmin disks, <2 4 can be 
calculated analytically for both theory classes. I find 
Q 4 = (8/15) 1 / 2 ' « 0.730 for the two MG theories, and 
Q 4 = (371-/16) 1 / 2 « 0.767 for all MI theories. These are 
a little larger [as expected since V(r) is increasing], and 
nearer each other than the Q values [as expected since 
the predicted V(r) of all theories merge asymptotically]. 
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